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Abstract 

We present results on the dynamics of the distorted diamond chain, S=l/2 dimers alternating 
with single spins 1 /2 and exchange couplings J\ and J3 in between. The dynamics in the spin fluid 
(SF) and tetramer-dimer (TD) phases is investigated numerically by exact diagonalization for up 
to 24 spins. Representative excitation spectra are presented, both for zero magnetic field and in 
the 1/3 plateau phase and the relevant parameters are determined across the phase diagram. The 
behaviour across the SF-TD phase transition line is discussed for the specific heat and for excitation 
spectra. The relevance of the distorted diamond chain model for the material Cus(C03)2(OH)2 
(azurite) is discussed with particular emphasis on inelastic neutron scattering experiments, a recent 
suggestion of one possibly ferromagnetic coupling constant is not confirmed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40. Gb, 78.70.Nx 
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I. INTRODUCTION 



The distorted diamond chain (DDC) is a one-dimensional (ID) quantum spin model with 
structure as shown in Fig. (Ha) and hamiltonian 

N/3 

H = | J2 S3 n+ iS3 n+ 2 + J\ (S3„S3 n+ i + S3 n+ 2S3„ + 3) + J3 (S3 n S3 n+ 2 + S3 n+ iS3„+3) | . (1) 

n=l 

This model with spins 1/2 and all couplings antiferromagnetic may be strongly frustrated 
owing to the triangular building blocks and has receiced increasing interest in the last decade 
for a number of reasons^: It has a rich quantum phase diagram as shown in Fig. QJb) (taken 
from Ref. 3)). Here and in the following we choose a representation with J2 = 1 as energy unit 
and Ji, J3 as variables. This representation emphasizes the symmetry of the model under 
exchange of J\ and J3. Three quantum phases have been discussed for the ground state of 
the model in zero magnetic field: For J 1; J 3 1 the ground state develops from the state 
with dimers in their singlet state on J2 bonds and nearly free spins between these dimers. 
The low energy sector is governed by an effective antiferromagnetic Heisenberg chain with 
N/3 sites resulting from the residual coupling between these spins and denoted as J c s in 
the following. This leads to the formation of a spin fluid (SF) phase with additional high 
energy excitations. For intermediate J\, J3 the ground state dimerizes, forming a twofold 
degenerate sequence of alternating tetramers and dimers (TD phase). Finally, for both J l5 J 3 , 
sufficiently large, the ground state is ferrimagnetic with e.g. a If J, structure of the unit cell 
of three spins (which satisfies J\ and J3 bonds and frustrates J2). These three phases can be 
clearly identified already in the symmetric model with Ji = J 3 in the regimes Ji = J3 < 1/2 
(SF phase), 1/2 < J x = J 3 < J m (TD phase) and J m < Ji = J 3 (ferrimagnetic phase) with 
J m « l.lOi. 

The generalization to the distorted diamond chain J\ ^ J 3 leads to an even richer be- 
haviour including e.g. the trimer Heisenberg antiferromagnetic chain on the J3 = axis 
with the standard HAF for J\ — 1 as limiting case. The DDC model can be seen as gen- 
eralization of the HAF with NN and NNN interactions: The critical point of this model at 
Jnnn/^nn ~ 0.2411 which marks the transition from a spin fluid to the dimerized phase 
is extended into the line between the SF and the TD phases in the DDC model and the 
Majumdar Ghosh point at Jnnn/^nn = 1/2 is extended into the line 1/2 < Ji = J 3 < J m 
with simple and exactly known dimerized ground states. In particular, the transition line 
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FIG. 1: (Color online) Structure (a) and phase diagram (b) of the distorted diamond chain model 
((b) from RefQ) 

between the SF and TD phases is a line of Kosterlitz-Thouless phase transitions. The point 
J\ = J3 = 1/2 is a point of particular high degeneracy. However, compared to the simple 
HAF chain with nearest and next nearest exchange, the DDC model has the dimer sub- 
system as an additional degree of freedom which dominates the high energy regime, but 
also, in the appropriate parameter regime, interacts with the low energy spin part. The 
ferrimagnetic part of the quantum phase diagram is another consequence of the combined 
influence of these two degrees of freedom. 

An external magnetic field polarizes the quasi-free spins and at a critical field H c \ produces 
a magnetization plateau at 1/3 of the total saturation magnetization. The plateau state 
corresponds to a fully saturated subsystem of spins 1/2 and all dimer s (J2 bonds) in their 
singlet state. Further increase of the moment requires breaking up at least one dimer with 
its large energy scale J2 = 1 and therefore a correspondingly higher field H > H c2 (end of 
the plateau). Finally complete saturation is obtained at the field H sat , given by 

^ = ~ + \{Ji + Js) + ~ {2(Jx - J3) 2 + \(Ji + -h? + 1 - ( h + Ja)} 2 • (2) 

The critical field H c i (beginning of the 1/3 plateau) is determined by the level crossing 
between the saturated state of the effective HAF and its 'ferromagnon' excitation with one 
unit of magnetization less. If the mapping to the effective HAF applies, this gives the relation 



H c i — 2J c ff. 



(3) 



Apart from the theoretical interest in investigating a model which allows to follow the 
variation between different quantum phases, the DDC model is of interest since it seems 



to describe reasonably well the compound azurite, Cu3(C03)2(OH)2. Azurite has been 
investigated in detail by static measurements (magnetization, susceptilibity, specific heat4) 
as well as by high field ESR- and the existence of the 1/3 magnetization plateau has been 
clearly established. From these experiments, this compound appears to be in the SF phase 
close to the phase transition to the TD phase. Recently, however, the possibility of one of the 
couplings Ji, J3 being ferromagnetic has been suggested from susceptibility and specific heat 
data£. Beyond the static properties investigated so far, the dynamics of the DDC and the 
material azurite in particular remain as a challenge to be understood both experimentally 
and theoretically: the characteristic feature of the model, namely the presence of two degrees 
of freedom with different energy scales and their mutual influence will show up most clearly in 
the energy spectra of the model. These are best investigated by inelastic neutron scattering 
(INS) experiments as clearly seen in recent work~. For a more complete description both 
of the DDC in the full phase diagram and of the results of INS experiments on azurite we 
present in the following results on the dynamics in the SF phase (section [TTJ) and in the TD 
phase (section [TTTI) . 

A perturbative approach can be applied to obtain results in the regime J±, J3 <C 1 (see 
Ref. |6|) as well as close to some special points in the phase diagram. Generally, however, 
for a quantitative description numerical calculations are required. We will present in the 
following the results of exact diagonalization, using both the Lanczos algorithm for systems 
with 24 spins and complete diagonalization (all eigenvalues) for 12 and 18 spins. The latter 
are necessary since the Lanczos algorithm gives only a limited number of the lowest energy 
levels in the subspace considered (in our case: S^ ot and wave vector k) which is not sufficient 
to cover the excited dimer subspace with its higher energies. Since the elementary cell has 
3 spins, our system sizes are restricted to 4, 6 and 8 elementary cells. It turns out, however, 
that this for many aspects is sufficient to obtain reliable results for the infinite system when 
a finite size analysis is carried through. 

II. DYNAMICS IN THE SPIN FLUID PHASE 

It is helpful to start the discussion from two well known limiting cases: 
(i) For J\ = J3 = the system reduces to A/3 independent dimers and A/3 free spins. In 
the ground state all dimers are in their sing let state and a 2 N ' Z fold de generacy due to the 
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free spins results. A magnetic field immediately saturates the free spin system leading to a 
magnetization of M sat /3 which remains constant until at a field H C 2 = J2 the dimers change 
to their triplet states saturating the system. This behavior to a certain extent remains valid 
on the symmetry line J\ = J3 where the total spins on all J 2 — bonds are independently 
conserved: the ground state as well as the 1/3 magnetization plateau in low field remain 
unchanged whereas the transition to full saturation is determined by the effective interaction 
which develops between two neighboring dimers in their triplet state and finally leads to an 
effective S — 1 chain. Qualitatively, for a large range of parameters the distorted diamond 
chain can be divided into two subsystems with clearly different energy scales, a low energy 
part of N/3 spins 1/2 and a high energy part of N/3 dimers. 

For small deviations from the independent free spin limit, J-y 7^ J 3 , the spin 1/2 subsystem 
develops some coupling by polarizing intermediate dimers and the 2 N ^ 3 fold degeneracy is 
lifted in favor of an effective Heisenberg chain with exchange J e g. In this regime, excitations 
of the DDC remain well separated: they are in the low energy regime with energy scale J c g 
forming the spinon continuum of the HAF with N/3 spins in the Brillouin zone of the full 
DDC (lattice constant a, reciprocal lattice vector r = lnj a) or in the high energy regime 
with energy scale J2 = 1 corresponding to the excitation of a dimer to its triplet state and 
developping into a dispersive band with width Adimer due to the coupling to the low energy 
spin subsystem. We will not consider in the following states with more than one excited 
dimer. 

(ii) For J\ = 1, J3 = the system reduces to the Heisenberg antiferromagnet with N spins, 
forming a spinon continuum in the Brillouin cell with reciprocal lattice vector 3r, energy 
scale 1 and no additional high energy excitations. For < J\ < 1 we have a trimerized 
Heisenberg chain and the spectrum is obtained by folding back the spinon continuum to the 
smaller Brillouin zone corresponding to lattice constant a. This results in three excitation 
branches (actually continua) which fill the energy range up to e = ir with small (for J\ 
slightly less then 1) gaps between them and an alternating sequence of minimum, maximum 
and minimum at wavevector k = n (in the following we will use exclusively the Brillouin 
zone with reciprocal lattice vector r, corresponding to the full DDC). With increasing 1 — J±, 
these trimer bands develop increasingly larger gaps, finally the continuum of the effective 
HAF emerges from the lowest band and the two upper bands conspire to give the dimer 
excitations decorated by continua of low energy spinon excitations. 
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Using this frame the lowest excitations of interest in the following are easily described: 

(i) the spinon continuum of the effective chain, 

(ii) the band with one excited dimer above the spinon continuum, 

(hi) the 'inverted ferromagnon' i.e. the saturated effective HAF with one spin deviation 
(3^ = ^/3-1), and 

(iv) the band with one excited dimer above the saturated effective HAF (one dimer in its 
triplet state (S£ ot = \N/Z + 1). 

The dispersion of excitations (ii)-(iv) is determined by hopping processes (spin deviations 
resp. dimer triplets moving to neighboring sites due to the residual interactions). To lowest 
order these processes result in a cosine dispersion and we introduce as notation for (ii) 

ir(*) = 1 + 4L + I A§L cos k. (4) 

More precisely, this excitation is not a single band but a continuum due to the spinon 
continuum of initial states; however, we will only be able to discuss the lower edge of this 
excitation and therefore simplify the notation using Eq.(j3]). Excitations (iii) and (iv) are the 
relevant excitations above the 1/3 plateau, we therefore use a notation giving their energies 
in finite magnetic fields relative to the plateau ground state with £f ot = |iV/3: 

eferrom(k) = ~Z A fcrrom (l + COS k) + H - H cl (5) 

= 1 + \ (Ji + J*) + S£L + \ A dtL(l + cos A;) - (H - H cl ). (6) 

The quantities A and 5 give the widths, resp. the nontrivial contributions to the minimum 
energy (at k — 7r) of the corresponding bands. The cosine dispersion of course is only valid 
in lowest order and will change to a more complicated expression for real systems. 

In the model of an effective HAF for the low energy regime its exchange constant J c g 
determines the low energy spinon (i) and the inverted ferromagnon spectrum: Af crrom = 2 J e g. 
Combined with e^^. it is also sufficient to give the range of the plateau phase and to 
characterize its dynamics: In the presence of a finite field, spectra are identical to those 
without field except for the shifts and splittings due to Zeeman energies. This establishes 
states with an increasingly larger total spin 5t Q t (i n their maximum z— component) as ground 
states. The plateau begins at the field H c i when the S tot = \N/?> level (saturated state of the 
quasi-free spin subsystem) is forced below the lowest 5t t = — 1 level ('ferromagnon' 

band at wavevector n) by the external field, leading to H cl = 2J cff . The lowest excitation 
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for the plateau dynamics close to the field H c \ then is the ferromagnon of Eq.(jSJ). When the 
field is increased across the plateau regime, the Stot = \N/3 + 1 level (one excited dimer on 
top of the saturated quasi-free spin subsystem) lowers its energy, crosses the ferromagnon 
excitation and finally is responsible for the end of the plateau at the upper plateau field H C 2 
implying 

H c2 = H cl + l+ l -{J l + J,) + 5^ CT . (7) 

The parameters determining the spectra can be calculated in perturbation theory in J 1; J 3 
and to lowest order are determined by the level spectrum of the general (J\ 7^ J 3 ) tetramer 
with 4 spins ... 3. This spectrum includes the lowest order information about J e gf in the 
singlet-triplet splitting of spins and 3 and about in the amplitude for the process 

I s) — > I J, t + ) (s and t + are noninteracting dimer states) which determines the propagation 
of an excited dimer triplet. The results to lowest order in Ji, J 3 are: 

2J e ff — ^ferrom = 2A ^ n pr = (Ji — J 3 ) 2 

^dimer = _ ^) 2 

H c2 = 1 + i(Ji + J 3 ). (8) 

J e fj has been calculated in straightforward perturbation theory up to fifth order-, based 
on the splitting of the general tetramer into singlet and triplet states we have obtained 
a result which accounts partly also for higher orders and allows reasonable estimates for 
\Ji — <h\ ^ 1, but finite J\ + J 3 : 

^cff = \{{i-h + J 3 ~ I) 2 + 3 {J\ - J 3 ) 2 }^ + h + J 3 - {1 + (Ji - J 3 ) 2 }~ 2 }- (9) 

In lowest order perturbation theory the relevant parameter, in addition to the energy scale 
set by J 2 and to J\ + J 3 , is the exchange of the effective HAF determined by ( J\ — J 3 ) 2 and 
many characteristic quantities of the DDC would be related by simple numerical factors if the 
mapping were perfect. Whereas these perturbational results allow to discuss the dynamics in 
principle, Ji, J 3 values of interest for the bulk of the phase diagram as well as for a material 
such as azurite are beyond the validity of perturbation theory. We therefore present in the 
following results from the numerical approaches described above. This will allow us to follow 
the essential aspects of the dynamics in the intermediate regime, i.e. through all of the SF 
phase. In Figs. [2] and [3] we show excitation spectra for three sets of exchange parameters: 
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set (a) represents the case of small couplings, set (b) is for a point in the phase diagram 
close to the SF to TD transition thought to be qualitatively representative for azurite^ and 
set (c) shows results for the case of one coupling ferromagnetic. The data are obtained by 
diagonalizing chains with 24 spins using the Lanczos algorithm. 

The low energy excitation spectra in zero external field are shown in Fig. [2J Here and 
in the following, energies are in units of J 2 and wavevectors in units of 7r/4 (for N = 24, 
resp. 7r/3 for N = 18). The spectra include all levels with S to t = 0, 1 of the spinon type in 
the low energy subsystem of 8 spins and in addition the lowest excitations for S tot = 2 (for 
completeness) and Stot = 3. The latter band of excitations is the 'inverted ferromagnon' and 
a cosine dispersion approximating the data points is shown for qualitative comparison to 
the effective model with its exact cosine dispersion. Fig. [3] shows the excitation bands in a 
magnetic field H c \ (beginning of the 1/3 plateau). In magnetic field two Zeeman components 
of the 5t t = 4, k = level are relevant: The Sf ot = 4 component turns into the plateau 
ground state, whereas the S£ ot = 3 component becomes the top of the inverted ferromagnon 
band, it is identical to that of Fig. [2] (apart from Zeeman shift) and now the lowest excitation. 
In addition, Fig. [3] shows the lowest excitation band with Stot = 5 which requires breaking 
one dimer (J 2 ) bond. Cosine dispersions as approximation to the data points are included 
for these two bands. For completeness we also show the first full band with Stot = 4 above 
the plateau ground state. 

Among the data shown, the excitations of interest from an experimental point of view 
(with large transition matrix elements for e.g. INS) are the spinon continuum in zero field 
and the inverted ferromagnon band as well as the excited dimer band in the plateau field. 
In addition, in zero field there will be transitions with energy of the order of J2 to an excited 
dimer band with Stot = 1 resulting from breaking a J2 bond on top of the effective chain 
groundstate. This is more difficult to deal with than the excited dimer excitation shown in 
Fig. [3] which is on top of the less complex saturated effective chain state. We will discuss 
these excitations below, based on calculations of all eigenvalues of a N = 18 chain (see table 
III). States with Stot > 1 in zero field as well as the band with Stot = 4 above the plateau 
ground state will be only weakly excited in INS and analogous experiments: in particular 
states in the S tot = 4 band are obtained from the saturated state by a virtual excitation 
I s T) ~^ \t+ 0- They have an excited dimer and an overturned spin (compared to the 
saturated state) in the low energy subsystem and thus require two spin flips to be excited. 
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FIG. 2: (Color online) Low energy spectrum of the DDC for N = 24, H = 0, J2 
Ji = 0.3, J 3 = 0, (b) Ji = 0.7, J 3 = 0.3, (c) Ji = -1.0, J 3 = 0.4. 
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FIG. 3: (Color online) Excitation spectrum of the DDC above the 1/3 plateau for N = 24, H = H c i, 
J 2 = 1.0 and (a) J x = 0.3, J 3 = 0, (b) J 1 = 0.7, J 3 = 0.3, (c) Ji = -1.0, J 3 = 0.4. 
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We now discuss how the dynamics changes with varying exchange constants: Set (a) 
shows the behaviour typical for the weakly coupled DDC: the bands are well separated in 
energy and the cosine dispersion is nearly perfect. Set (b) displays what is expected for a 
material such as azurite: in zero field a spinon continuum should be clearly visible whereas 
in the plateau regime two separate bands dominate the picture. The cosine approximation 
to the dispersion is less applicable, actually the spectrum of the inverted ferromagnon is 
close to linear for smaller wave vectors. For the set (c) which serves as an example for the 
alternative suggesting one ferromagnetic coupling^, the dynamics in zero field is seen to be 
surprisingly close to that of set (a). This may explain the emergence of the ferromagnetic 
alternative from a discussion of static quantities. However, these two sets lead to strongly 
differing dynamics in finite field as seen by comparing Fig. [3] (b) and (c): The standard 
antiferromagnetc model (b) implies two well separated bands with rather small widths, 
whereas the partly ferromagnetic alternative (c) is characterized by an overlap of the two 
bands and a strong dispersion of the excited dimer band. 

For a semiquantitative discussion of the low energy dynamics of the DDC the mapping 
to the model of an effective HAF is rather useful. Therefore we discuss shortly the quality 
of this mapping for H = in Fig. H]for the two parameter sets (a) J\ = 0.3, J3 = and (b) 
J\ = 0.7, J3 = 0.3. We compare the numerical spectrum for N = 24 (dots) with the energies 
of the N = 8 HAF with unity exchange constant (open circles). The energies of the DDC 
have been scaled by an effective exchange constant J e g, chosen to reproduce the (N = 8) 
maximum spinon energy at k = |, S tot = 1. In (b) energy levels of the N = 8 HAF which 
are beyond the range of the Lanczos calculation for the DDC model have been omitted from 
the plot to obtain a clearer picture. Whereas the mapping for the small parameter values in 
set (a) is nearly perfect throughout all of the spectrum, substantial deviations are seen for 
the parameter set (b): the low energy spinon part is still reproduced well by the effective 
model, but the high energy part, in particular the ferromagnon band with Stat = 3 is very 
different both in energy and in dispersion. A cosine dispersion is only a rough approximation 
to the spectrum. 

Thus a quantitative experimental investigation of the dynamics will contribute substan- 
tially to locating the position of a specific compound in the phase diagram of Fig. [TJ(b). For 
a quantitative overview (and possibly use in determination of coupling parameters from ex- 
periment) we reduce in the following the information in these spectra to a few characteristic 
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FIG. 4: (Color online) Spectrum of the effective HAF (N=8) vs spectrum of the DDC (N=24) for 
(a) Ji = 0.3, J 3 = and (b) J x = 0.7, J 3 = 0.3. 

numbers to be presented in tables H and [III below. From the numerical data we have calcu- 
lated values for the quantities determined by the effective exchange between the quasi-free 
spins. We give in table [J numbers for J e g determined from the maximum spinon energy at 
k = 7r/2, 5tot = 1 (when multiplied by 1.7964.., the corresponding energy in the N=8 HAF 
chain, these numbers lead back to the energy for the DDC model) and for the S tot = 1 spinon 
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at k = 7r (gapped due to discreteness). Further we give the width of the ferromagnon band 
(which would be 2J e fj if the mapping to the effective model were perfect) and the width of 
the dimer band above the plateau with S tot = N/6 + 1. In table [Til we give numerical values 
for the characteristic magnetic fields, i.e. beginning (H c i) and end (-£^2) of the plateau as 
well as the saturation field H sat . In the standard case (actually some exceptions exist close 
to the phase transition line) H c \ is identical to the ferromagnon width from table [TJ In 
table [TTJ we also give the energy scale which is relevant for an application of the numerical 
results to azurite: Using the experimental number if sat = 33T the value of the coupling J 2 
is calculated from Eq. (j2J) and the values given in table [TTJ (in both meV and T) may serve 
to obtain energies and fields applying to azurite in standard units. In tables ITI and ITTI three 
regimes of the SF phase are covered: (a) values along the J3 = axis (i.e. for Heisenberg 
trimer model) (b) values along a diagonal path which appears as the most interesting one for 
discussing azurite and (c) tow examples for ferromagnetic coupling. We will discuss below 
the possibility of such an interaction from the point of view of inelastic excitations. Values 
along the line J\ = 0.6, passing through the phase transition at J3 ~ 0.364, will be presented 
in section ITVl 

If the mapping to an effective low energy Heisenberg chain were perfect, it would imply 
relations between three different quantities which are all determined by J e s: 

(i) the maximum of the effective one spinon dispersion at wave vector tt/2: e(ix/2) = 7r/2J e fj, 

(ii) the critical field which determines the beginning of the plateau: H c \ = 2 J e g, 

(iii) the width of the effective 'inverted ferromagnon' as the lowest excitation at the be- 
ginning of the plateau (in fact, its minimum defines H c i): €f m (0) — ef m (7r) = 2J c g. This 
is automatically equal to (ii) from the definition of H&, but the cosine dispersion is an 
additional independent property. Since the mapping is only approximate, these quantities 
differ as is seen in the numerical data and the differences characterize the quality of the 
mapping. Actually there are more possibilities to extract J e fr from the numerical data such 
as the energy of the lowest spinon singlet at k = it and the ground state energy (suit- 
ably extracted from the energy of the saturated subsystem state), but numbers from these 
approaches essentially confirm the picture as it has emerged from the tables above. The 
essential conclusion for the real DDC is that the effective coupling J e g depends on energy. 

The quantities in tables [TJ and [TT] have been calculated for N = 24, but a comparison with 
results for N = 12 and N — 18 shows that finite size effects are very small, indeed negligible 
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TABLE I: DDC parameters related to the effective interaction (see text) 



couplings Ji , J3 




H C 2 


H sa x, 


energy scale J2 in 
meV resp. T 


0.02, 0.00 


4.12 10~ 3 


1.0100 


1.0102 


4.16 meV = 32.7 T 


0.30, 0.00 


0.121 


1.143 


1.200 


3.50 meV = 27.5 T 


0.60, 0.00 


0.480 


1.245 


1.500 


2.80 meV = 22.0 T 


0.80, 0.00 


0.795 


1.272 


1.740 


2.41 meV = 19.0 T 


0.70, 0.30 


0.529 


1.390 


1.627 


2.58 meV = 20.3 T 


0.65, 0.25 


0.452 


1.371 


1.569 


2.68 meV = 21.0 T 


0.60, 0.20 


0.388 


1.342 


1.512 


2.78 meV = 21.8 T 


0.40, 0.00 


0.220 


1.183 


1.29 


3.26 meV = 25.6 T 


0.40, -1.00 


0.478 


0.861 


1.234 


3.40 meV = 26.7 T 


0.02, -0.40 


0.091 


0.824 


0.868 


4.84 meV = 38.0 T 



TABLE II: DDC parameters related to the characteristic fields and unit of energy for azurite (see 
text) 
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excited dimer dispersion for H = H_cl 
combine N=18 and N=24: S A z_tot = N/6 +1 



A 



12 3 4 

wave vector 

FIG. 5: (Color online) Excited dimer band of the DDC (N=24 and N=18) above the 1/3 plateau 
ground state (in a magnetic field H c {) for J\ = 0.7, J3 = 0.3 and J\ = 0.8, J3 = 0.2 

to the accuracy given. This is due to the fact that e.g. J e s is determined by comparing two 
finite systems, the HAF chain with N = 8 and the DDC with N = 24. This is illustrated in 
Fig. [5] where the excited dimer bands for two different sets of couplings are shown combining 
results for N = 24 and N = 18 in the same graph. We conclude that the the effective 
parameters considered so far can be reliably determined for the infinite chain. The situation 
is different for the last quantity of interest, the dynamics of the one dimer excitations for 
H = 0. These excitations originate by exciting one dimer from singlet to triplet on top of the 
spinon continuum of the effective HAF chain (instead of on top of the saturated effective 
HAF chain for Adimer as given above). In the following we give results on this zero field 
branch as far as numerically possible: 

At zero field, the excited dimer states form a continuum as well and the infinite chain 
limit can only be obtained by extrapolation in 1/N. With iV = 12 and N = 18 as the 
only available numbers of spins the extrapolation can only be done for wavevectors k — 
and k = it. The lowest energy levels in our finite systems result from coupling of spinons 
at wavevectors (singlet only) and it (singlet and triplet) to the dimer triplet, i.e. we 
get a band of one singlet, three triplets and one quintuplet. However, the separation both 
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couplings J\ , J3 


e dimer (fc = 0) 


e dimcv (k = vr) 


0.30, 0.00 


0.861 


0.893 


0.60, 0.05 


0.988 


0.978 


0.60, 0.25 


0.720 


0.711 


0.65, 0.25 


0.854 


0.793 


0.70, 0.30 


0.864 


0.776 


0.40, -1.00 


1.226 


1.255 



TABLE III: Energies e(k = 0) and e(k = ir) in the DDC lowest excited dimer band (onset of the 
continuum). Values extrapolated from TV = 12 and N = 18 to infinite N 

between these multiplets and to the higher levels is due to the finite size, whereas in the 
thermodynamic limit a continuum with energy of the order of J 2 will result. This is to some 
extent reflected in the results of finite size extrapolation of the multiplet energies related 
to one excited dimer: There is a clear tendency for the energies to converge to the same 
value. Thus only one energy in this high energy subspace can be given reliably, no reliable 
information can be obtained in this approach about splitting into bands. In table II III we 
give the energy of the lowest state in the continuum of excited dimers obtained this way for 
a number of coupling constants. 

III. DYNAMICS IN THE DIMER-TETRAMER PHASE 

The ground state in the tetramer-dimer phase is twofold degenerate and develops from 
the ground state on the symmetry line J\ = J3. On this line the two ground states can be 
written down explicitly (even for finite systems with an arbitrary even number of cells of 3 
spins). They are given by the alternating sequence of the dimer singlet S and the lowest 
tetramersinglet T. This allows the two equivalent configurations 

...S T S T S T... and ...T S T S T S... (10) 

describing the two degenerate ground states. 

The lowest excitations above these ground states are obtained as solitons which are defined 
by gluing together the two degenerate ground states in a localized region on the chain. This 
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gives e.g. the state 

...S T S * S T... (11) 

where * denotes a free spin. A soliton is possible only with two dimer singlets adjacent 
to each other and a free spin between them whereas a configuration . . . T T . . . obviously 
does not exist. On the symmetry line J\ = J3 there are N/3 degenerate localized one- 
soliton configurations. They start to propagate and to form a soliton band for J\ ^ J3. The 
properties of a single soliton can suitably be investigated for chains with an odd number 
of cells N/3 when periodic boundary conditions require the existence of one soliton in the 
ground state. The hopping process 

T S * -> * S T (12) 

with amplitude t leads to the propagation of solitons and the formation of a ground state 
band with energy 

e sol (A;) = E + 2tcos2k. (13) 

Thus, from numerical calculations for an odd number of cells the hopping amplitude 
is easily determined even when only small systems are available. For infinite chains with 
periodic boundary conditions and an even number of cells iV/3 which possess two degenerate 
ground states and enforce an even number of solitons, the low energy spectrum is dominated 
by the two soliton continuum emerging from independent propagation of two solitons with 
wavevectors \k + k\ and \k — k±, i.e. 

e 2so \k) = 2E + At cosk cos2A;i. (14) 

Numerically obtained soliton spectra are shown in Fig. [6] for (a) N = 15 and (b) N = 18 
for the point J\ = 0.6, J3 = 0.55 close to the symmetry line in the phase diagram. e sol (fc) in 
Fig. E^a) clearly shows the cos 2 A; dispersion, whereas the dispersion of e 2sol (fc) in Fig. [6](b) is 
somewhat more complicated due to the small size of the system. Also shown are the three, 
resp. four-soliton bands demonstrating the clear division of the spectra into distinct solitons 
bands for this nearly symmetric set of couplings (the zero of energy in Fig. EJ^a) is taken 
from the noninteracting limit). The situation is analogous to the Ising chain with small 
transverse interactions, the system where the dynamics of magnetic solitons was discussed 
firstr^i. Slightly different, the soliton spin 1/2 here is a real spin 1/2 which can be attributed 
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(a) J_2 = 1.0, J_l = 0.6, J_3 = 0.55 : solitons in TD phase (N = 15) 
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(b) J_2 = 1.0, J_l = 0.6, J_3 = 0.55 : solitons in TD phase (N = 18) 




wave vector 



FIG. 6: (Color online) Soliton spectrum of the DDC ((a) N=15, (b) N=18) in the TD phase for 
J\ = 0.6, J3 = 0.55. The N=15 spectrum is shown in the complete Brillouin zone 0...27T to 
demonstrate the cos 2k dispersion of the single soliton. 

to the free electron of the Cu 2+ ion between the two dimers forming the domain wall. For 
small numbers N/3 the soliton spectrum clearly shows the effects of the different symmetry 
of singlet and triplet. 

From the one soliton data for N = 15 and J\ = 0.6, J3 = 0.55 the hopping amplitude is 
deduced as t ~ 0.0076. For the two soliton data for N = 18 the corresponding calculation 
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has to include the possibility of two neighboring solitons as well as the resulting symmetry 
effects and gives a somewhat higher value, t ~ 0.0097. The deformation of the zero order 
wave function due to neighboring solitons is strongest for small systems which explains the 
difference. However, for N = 24 (when only the lowest 2 soliton band is accessible in Lanczos 
calculations) we obtain t w 0.0076 ± 0.0002, identical to the one soliton result for N = 15 
within the uncertainty resulting from matching the cosine dispersion for the different wave 
vectors. 



IV. CROSSING THE PHASE TRANSITION LINE 



In this section we present some data for the specific heat and for the spectrum of low- 
lying excitations in order to approach the behavior of the system when its coupling constants 
change between well defined end points, one in the SF phase, the other one in the TD phase, 
thus crossing the phase transition line. Evidently, owing to the small system sizes accessible 
only in our calculations, we cannot claim that these data describe correctly the most inter- 
esting aspect, namely the critical behavior; on the other hand our data for both the specific 
heat and the spectrum of low-lying excitations set a reasonable frame for the transition 
regime, to be filled by more detailed calculations later. In the following we present results 
for the DDC on the line J<i = 0.6 for varying J 3 . As discussed above, the phase transition 
along this line is of Kosterlitz Thouless type and can be considered as a generalization of 
the phase transition in the HAF with both nearest and next nearest neighbor exchnage. We 



therefore have applied the procedure of Ref. 
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to determine the critical coupling and find 



that the phase transition occurs at J3 ~ 0.364. 

Fig. [7J shows a sequences of specific heat data varying J 3 for fixed J± for three different 
values of the magnetic field. The data are obtained from the full spectrum for the N = 18 
chain and therefore cover reliably the complete temperature regime although critical prop- 
erties near the critical coupling J 3 « 0.364 will appear smeared out. In all diagrams we 
use a logarithmic temperature scale adequate to the strongly different energy scales. For 
all magnetic fields the specific heat exhibits a high temperature peak at T « 0.5, whereas 
the low temperature properties reflect the structure of the system: For H = and low tem- 
perature (Fig. [Th.) the SF phase is characterized by a continuously increasing contribution 
to the specific heat, the remnant of the effective Luttinger liquid. With increasing J3 this 
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(a) Specific heat for H = 0, J_3 = 0.6 
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(c) Specific heat for H = 1.3, J_l = 0.6 
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FIG. 7: (Color online) Specific heat of the DDC (from all levels for N = 18) for H = (a), H = 0.6 
(b) and H = 1.3 (c), fixed J\ = 0.6 and varying J3 through the phase transition. 
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FIG. 8: (Color online) Low-lying spectra (5 to t = (red/full) and 5 to t = 1 (blue/open)) of the 
DDC (N=24) for J x = 0.6 and (from top to bottom) J 3 = 0.55, 0.44, 0.32, 0.05. The energy range 
for each spectrum is 0.5. 

contribution develops gradually into the gapped contribution of the TD phase with charac- 
teristic shoulders on both sides of the phase transition at J3 ps 0.364. For H = 0.6 (Fig. [7b) 
the DDC always is in the gapped plateau regime and the specific heat shows little variation 
with the coupling. For H — 1.3 (Fig. [7b) the DDC is close to saturation and the differences 
between the two phases, resulting from the energy spectra above the plateau gap, become 
apparent again. In particular, the well defined grouping into soliton bands at J3 = 0.55 
leads to a strong low temperature peak which actulaly develops continuously from lower J 3 
values. It would be interesting to see, using more powerful numerical methods, whether the 
development of this peak in the nearly saturated case shows critical properties. 

Fig. [H] shows a sequence of spectra of low-lying excitations (S tot = 0, 1) in zero magnetic 
field obtained by the same procedure as the spectra shown in sections [III and IIIII Qualita- 
tively, the transition from the SF phase with its gapless spinon continuum to the gapped 
soliton spectra in the TD phase is clearly seen qualitatively, as far as possible for the limited 
size of the system: with increasing J3 the second degenerate ground state at k = ir emerges, 
the spinon continuum is compressed into the soliton band and the gapless character dis- 
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appears. Excitation spectra in finite magnetic field, in particular for fields in the plateau 
regime, on the other hand do not show specific variations but rather continuous changes 
across the phase diagram without prominent features close to J3 = 0.364, the location of 
the phase transition in zero field. This illustrates, as expected, that the signature of the 
phase transition in the dynamics is limited to zero external field and low-lying excitations 
wehreas the 1/3 plateau as well as the accompanying excitations are continuous across the 
phase diagram. 

V. CONCLUSIONS 

For the S = 1/2 distorted diamond chain in both the spin fluid (SF) and the tetramer- 
dimer (TD) phases we have calculated the spectra of low-lying excitations and the specific 
heat. We have used both full numerical diagonalization (for chains with up to 18 spins) and 
the Lanczos algorithm (for chains with up to 24 spins) and have discussed the results in 
relation to approximate analytic approaches. Except close to the SF-TD phase transition 
results for our small systems are shown to represent the thermodynamic limit. Our calcula- 
tions are for arbitrary value of the external magnetic field, results are mainly given for zero 
field and for fields corresponding to the 1/3 plateau regime. 

In the SF phase the low energy spectra can be related to a Heisenberg antiferromagnetic 
chain with effective interaction J c g. For parameters beyond the validity of a perturbative 
approach, this effective interaction has to be allowed to be energy dependent. The lowest 
excitations in the plateau regime are the inverted ferromagnon and the propagating sin- 
gle dimer triplet excitation with, however, partly strong modifications of the corresponding 
cosine dispersions. The values of the characteristic parameters (J e s, extent of the plateau 
regime, widths of the cosine bands) are given for typical paths crossing the SF phase. These 
data should allow to decide whether a material such as Cu3(C03)2(OH)2 (azurite) is suffi- 
ciently well described by the DDC model and, if so, to determine the corresponding cou- 
plings. The standard assumption for azurite is to take all couplings as antiferromagnetic 
and we have shown that the spectra of low-lying excitations exhibit large and characteristic 
changes when the possibility of one ferromagnetic coupling is introduced. We therefore ex- 
pect that our data will allow to interpret quantitatively experimentla data on azurite. This 
refers in particular to the results of inelastic neutron scattering experiments^. Considering 
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the present status of such investigations, our results do not confirm the conclusion of at 
least one ferromagnetic coupling in azurite. Generally, our results lead us to describe the 
following signatures when ferromagnetic couplings are present: 

(i) Whereas the dimer width is roughly 1/2 of the ferromagnon width for af couplings (as 
suggested by perturbation theory), for ferromagnetic couplings these widths tend to become 
equal. 

(ii) The sign of the couplings has a marked influence on the relative appearance of the ferro- 
magnon and the excited dimer band: For couplings Ji, J3 = (—1, 0.4) these bands above the 
1/3 plateau overlap, whereas for (—0.4,0.02) ferromagnon a swell as dimer width become 
very small and correspondingly if cl becomes much smaller than in perturbation theory. 

The low-lying excitations in the TD phase with its twofold degenerate ground state are 
shown to be solitons. The width of the one soliton band as determined from the N = 15 
chain not too far from the symmetry line J\ = J3 reproduces well the soliton bands in the 
N = 24 chain and therefore gives reliably the tunneling amplitude for the soliton propagation 
in the thermodynamic limit. 

We have also shown spectra as well as the specific heat on a line across the SF-TD phase 
transition. Although the small systems accessible to us do not allow to discuss critical 
properties of the DDC close to this Kosterlitz-Thouless transition, the variation of the 
dynamical properties through the transition become clear. In particular, only the low energy 
properties, determining the behavior of the system at zero field, carry the signature of the 
phase transition. 

Generally, the DDC has many features in common with the antiferromagnetic Heisenberg 
chain with nearest and next-nearest exchange and the SF-TD phase transition is of the same 
type as the KT transition at Jnnn = 0.2411...Jnn in this system. On the other hand, we 
have shown that the additional degrees of freedom, resulting from the possibility to excite 
the J2— dimers to the triplet state, show up clearly in the dynamics. We leave to the future 
to investigate the influence of these degrees of frededom on the phase transition using more 
powerful anlytical and numerical methods. 
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